CPT-2004/P.006 
ROMEl-1370/2004 
ROM2F/2004/02 



Topological susceptibility in full QCD 
with Ginsparg- Wilson fermions 

L. GiustP, G.C. Rossi^ M. Testa'^ 

°- Centre de Physique Theorique, CNRS Luminy, Case 907, F-13288 Marseille, France 

^ Dipartimento di Fisica, Universita di Roma "Tor Vergata" 
INFN, Sezione di Roma 2 
Via della Ricerca Scientifica, 1-00133 Roma, Italy 

^ Dipartimento di Fisica, Universita di Roma "Xa Sapienza" 
INFN, Sezione di Roma "La Sapienza" 
P.le A. Moro 2, 1-00185 Roma, Italy 



Abstract 

We show that, if the formula for the topological charge density operator suggested by 
fermions obeying the Ginsparg- Wilson relation is employed, it is possible to give a precise 
and unambiguous definition of the topological susceptibility in full QCD, XtL^'j for finite quark 
masses on the lattice. The lattice expression of xIl^ looks like the formal continuum one, in 
the sense that no power divergent subtractions are needed for its proper definition. As a con- 
sequence, the small mass behaviour of xIl^ leads directly to a multiplicative renormalizable 
definition of the chiral condensate that does not require any power divergent subtraction. 



1 Introduction 



In this paper we discuss the definition and the properties of the topological susceptibility in full 
QCD with massive quarks extending the results of ref. Using arguments based on anomalous 
flavour singlet Ward-Takahashi identities (WTI's), we prove that, if the formula of the topological 
charge density, Q, suggested by Ginsparg-Wilson (GW) fermions |2l E] is employed, the full 
QCD topological susceptibility 

xlf = / d^x{Qix)Q{0)) (1) 

does not need any power divergent subtraction at finite non-vanishing values of the quark (pion) 
mass. Furthermore it vanishes linearly in the quark mass with a coefficient which turns out to be 
(proportional to) the chiral condensate, as in the formal continuum limit. 

The interest of these results lies in the fact that one can exploit the absence of power divergent 
mixings in the continuum-looking lattice formula to extract the value of the chiral condensate 
with no need of performing any dangerous power divergent subtraction. 



2 Generalities on Ginsparg-Wilson fermions 

Regularizing the fermionic part of the QCD action using GW fermions offers the great advantage 
that global chiral transformations can be defined, which are an exact symmetry of the massless 
theory, as in the formal continuum theory. This is a consequence of the relation |2j 

j^D + = aDj^D , (2) 

where D is the Dirac operator and a is the lattice spacing. Eq. (j2I) implies the invariance of the 
massless fermion action under the transformations 

6{^fj = \%iP, = i^-i^Xf , / = 0,l,...,iV|-l, (3) 

where \^ = \ and the A-'^'s, / 7^ 0, are flavour matrices ^. In the flrst of the equations above we 
have introduced the definition 

75 = 75(1 -aD) (4) 

with the properties 

75 = 75 ' 75 = 1 ■ (5) 

Eqs. (jSJ may be interpreted as the lattice form of the continuum chiral transformations. The 
Neuberger operator |^ satisfies the GW relation, has the correct continuum limit and is local, 
though not ultra-local [S] . Another solution of the GW condition is provided by the fixed-point 
fermionic action of refs. [Z|IH]- 

^We use the normalization ir{\^ X') = <5-^»/2, [\^ ,\^] = if^^'^X'^, so that {\^ ,\<'} = d^^'^X'' + S^^l/Nf , f,g,h = 
l,...,Nj-l. 
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In a GW regularization the lattice QCD fermion action with Nf massive flavour can be written 
in the form ^ 



Sf 



where 



(fx J2 

r,s=l 



[X) \ (D6rs + P-Ml,P_ + P+MrsP+) V'l (X) 



Vl±75), P± = ^(l±75), 



(6) 



(7) 



2' ^ 2 

ip (ip) is an A^j-dimensional column (row) vector in flavour space and M = diag (mi, . . . , rriMf)- Sp 
is invariant under the UiiNf) x Ufi{Nf) global transformations 



■0L 



with Ul,r G U{Nf)L,R and 



^R 



P+^ 



IpL 
i^R 



i'R 



^rU] 



R ' 



■i^P-, 



(8) 



(9) 



if at the same time M UlMUj^. In the following we shall restrict ourselves to the flavour vector 
symmetric case nir = m, with r = 1, . . . ,Nf. As a consequence of the exact chiral invariance of 
the massless GW regularization, no additive quark mass renormalization is required. The action is 
0(a)-improved, since no chiral invariant operators of dimension d = 5 can be constructed. In this 



work we will consider the bilinear scalar and pseudo-scalar quark operators (/ = 0, 1, 



S^ix) =i){x)\f 
pf{x) = V;(x)A^75 



-In), 
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-1) 
(10) 

(11) 



The "rotation" (1 — |D) of the quark field ijj in the above equations leads to definitions of scalar 
and pseudo-scalar quark densities which have the correct chiral transformation properties, like in 
the formal continuum theory, and only need a (logarithmically divergent) multiplicative renormal- 
ization. Furthermore the operators pOj) and are automatically 0(a)-improved. 

In a GW regularization the gauge anomaly is recovered a la Fujikawa |31 The fermion 
integration measure is not invariant under Ua{^) transformations (eqs. (jH)) with / = 0), and the 
topological charge density 



a'^Qix) 



-Tr 
2 



75-D(x,x) 



(12) 



originating from the corresponding Jacobian, is related to the index of the lattice Dirac operator, 
D, by the equation [lITniin] 



n+ 



n_ 



index {D) = / d^x Q{x 



(13) 



^For short we use continuum looking notations with j d'^x replacing y^_^ 
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with n+ the number of zero modes with positive (negative) chirality ^. For a recent review 
on this subject see [TT] . 

3 The singlet Ward-Takahashi identities 

In the chiral hmit, the local anomalous flavour-singlet WTI's have the form 

d,{Al{x)d{y)) = 2Nf{Q{x)d{y)) - {6'/d{y)) , (14) 

where J^^{x) is the singlet axial current, O is any renormalized (multi-)local operator, concentrated 
at points y = = 1, . . . ,n} and 6^j^O is its local singlet axial variation. In eq. 1)14(1 we have 

not shown the exponentially suppressed terms coming from the fact that D is not ultra-local [Sj. 
They are of no relevance for the following arguments, as they vanish after integration. Assuming 
the absence of a Ua{^) massless Goldstone boson, the integrated form of the WTI's (|15() reads 



= 2iV^ y <fx {Q{x)0{y)) - {5\0{y)) . (15) 

Since the second term in the r.h.s. of eq. (jl5j) is finite, it follows that / d^x Q{x) is also finite, as it 
has finite insertions with any string of renormalized fundamental fields. Therefore Q{x) can only 
mix with operators of dimension < 4 and vanishing integral, hence only with d^J^^{x). No power- 
divergent subtractions with lower dimensional operators (such as the pseudo-scalar quark density) 
are required This is a very distinctive feature of GW fermions with respect to standard Wilson 
fermions which is directly related to the absence of an additive mass renormalization ^. One can 
define finite operators Q and jP^ by the equations ^ 



Q{x) = Q{x) - ^d,A%x) , (16) 

Alix) = {l-Z)Al{x), (17) 

where Z is the mixing coefficient between Q and d^j^A^. The renormalized singlet axial WTI's then 
become (again up to exponentially small terms) 

d,{Alix)diy)) = 2Nf{Q{x)d{y)) - {6'/d{y)) . (18) 

Outside of the chiral limit the integrated singlet axial WTI's read 

= 2mJ d^x {P^{x)d{y)) + 2Nf J d^x {Q{x)d{y)) - J d^x {6^/d{y)) . (19) 

The extra term present in eq. (|16|) being a total divergence does not contribute to the integrated 
WTI's ()19|) . If we replace 0{y) with the local operator Q{0), we get 

= 2mJ d^x (P°(x)Q(0)) + 2Nf J d^x (Q(x)Q(O)) , (20) 



■^For alternative lattice definitions of Q see the papers of ref. 

"'a discussion of the singlet WTI's for Wilson fermions can be found in ref. 
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and, similarly, by inserting the multiplicative renormalizable operator P''(0) 

= 2mJ (fx (P°(x)P°(0)) + 2Nf J (fx (Q(x)P°(0)) - 2(S°(0)) . (21) 
Putting together eqs. H2U() and (PT|) . we obtain 

^ / d'x {Q{x)Qm = / d'x {p\x)p'm - J^is'^m . (22) 



(2A^/)2 J ' ' {2Nfy 

In the next section we show that the full QCD topological susceptibility, xfi^^ defined above is not 
affected by power divergences. 

3.1 Absence of power divergences in XtL^ 

The proof of the absence of power divergences in Xt"^^ is based on the study of the short distance 
beahviour of the two terms in the r.h.s. of eq. H22|) at small m. 

The argument goes through a number of steps. First of all we observe that thanks to the chiral 
properties of GW fermions, only power divergences of the type rr? jo? can possibly be present in 
eq. ()22|1. The second observation is that the m? /a? divergences separately affecting the two terms 
in the r.h.s. actually cancel each other. This is the result of the exact compensation between the 
(quadratically) divergent term arising in m? J d^x {P^ (x) (0)) , due to the short distance behaviour 
of the integrand, and a similar divergent term appearing in (S^), when one power of the fermionic 
mass term (brought down from the action) is inserted together with S^. The compensation follows 
from the fact that the short distance (perturbative) behaviour of the two correlators {P^{x)P^{0)) 
and {S^{x)S^{0)) are (in the massless limit) equal (up to a minus sign due to the presence of two 
extra 75 matrices in {P^{x)P^{0))), leaving behind a finite, computable contribution. 

We now make explicit and precise the line of arguments sketched above. 
1. Each term in eq. H22|l is even under the (non-anomalous) spurionic symmetry ^3] 

(23) 

V' -^tp' = -V'75 

where TZ^ is an element of the chiral group. Since only the identity operator, which is even 
under TZ^^, can contribute a divergent term in the two terms in the r.h.s. of eq. (|22|) . we 
conclude that only mP/a^ power divergences can be present, as they are even under m —m. 
In other words chiral invariance forbids power divergences like mjc? and rr? ja. 



2. If we order the terms contributing to eq. (|22j) in powers of m, we get 

4iV?XtT = -4m(5°(0)) +4m2 /d^x [(5°(x)5°(0)) + (P°(x)P°(0))l + (24) 

^ m=0 J l J m=0 



+ 



0(m^/m^) with/without S^SB 
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where odd powers of m can be present in the expansion only if chiral symmetry is spon- 
taneously broken. For the purpose of studying the structure of power divergences, we only 
need to examine terms 0(m) and 0(m-^). Higher order terms are at most logarithmically 
divergent ^. 

3. The first term in the expansion (|24|) is finite. In fact, i) owing to the exact chiral symmetry of 
the massless GW fermionic action, there cannot be any mixing between the identity and the 
operator 5" (with a divergent coefficient), because they transform in the opposite way 
under TZ^^; ii) the quantity m{S^{0))\m=o is not logarithmically divergent, as a consequence 
of the non-singlet WTI's. 

4. The sum of the next two terms is finite. To prove this result it is convenient to consider the 
set of WTI's {f = l,...,Nj -1; g,h = 0,...,Nj -1) 

= Jd'^z J d^xdf,{Al{z)P3{x)S''{0)) = 
= 2mJ d^z J d^x {pf{z)PS{x)S''{0)) + 

d^x {d{p9{x)S''{0)) - [ d^x {P3{x)5{s'^{Q)) , (25) 



where represents the operation of taking the axial variation with flavour index /. 

5. Combining the above WTI's, one gets in the chiral limit the soft-pion theorem (no sum over 
repeated indices, d^^^ ^ 0) 



m=0 



?TI=0 



+ 



d^x [{S^{x)S\<d)) + (P°(x)P°(0)) 
= F^Nf[ I d^x{7Tf\TEiPHx)S\0) + Sfix)P'^iO))\0) 

where Te means Euclidean time-ordering ^ and we have used the definition 



From the O.P.E.'s 



P^(x)S°(0) ~ 5^(x)pO(0) ~ ^P^(O) , 

X'^ 

pfix)SS{0) ~ J2 dff^'^^P^{<d) , 

h 

one concludes that the integrals in the r.h.s. of eq. (|26() are indeed finite. 



(26) 

(27) 

(28) 
(29) 



^Whether the logarithmically divergent terms proportional to might be reabsorbed by renormalizing m is a 
question that can be decided by a perturbative calculation. 

^For simplicity we do not employ a different notation for the operators in the matrix elements appearing in the 
r.h.s. of eq. 1261 . 



4 Final considerations 



A number of interesting consequences follow from the formulae ()18() and ()22|). 

1. In the full theory m^, ^ and there is no massless particle that can couple to . So it is 
immediately seen that XtL^ vanishes in the chiral limit (m — > 0). 

2. A formula for the (quenched) rj mass jl61 117j can be obtained starting from the Fourier 
transform of the WTI (|18|) at zero quark mass, if one chooses O = Q. For completeness we 
recall here the standard argument which goes as follows. First we observe that the C/a(1) 
variation of Q is zero. Taking the Fourier transform of eq. (|18j) with O replaced by Q, one 
gets (in the chiral limit) 



i J d^xe-*P%(^)l(x)g(0)) = 2Nf J d^xe-*P"(g(x)g(0)) . (30) 

In the limit Nf/Nc — > 0, where the rj' mass vanishes, the l.h.s. of eq. (|30|) is dominated at 
small p by the rj' pole, leading to the expansion 

i /d4xe-^P%(^0(x)Q(0)) , = lim mlF^, o ^\ + Q(p') • (31) 

If, as indicated in the above formula, the limit p — > is taken after the limit Nf/Nc 0, one 
ends up with the relation 



2N 



f 



lim lim [ d'^xe-'P''{Q(x)Qm, (32) 

P^ONf/Nc^oJ 



Nf/Nc=0 p^ONf/Nc^O 



where standard counting arguments JH] ensures that the l.h.s. of eq. ()32() has a finite limit 
as Nf/Nc 0. At this point if it is assumed that, taking the limit Nf/Nc ^ in the r.h.s. 
of ()31|) is equivalent to drop the fermion determinant, one arrives at the famous WV formula 



2Nf 



Nf/Nc=0 



(fx (g(x)Q(O)) . (33) 



YM 



Notice that in the limit Nf/Nc 0, the mixing coefficient, Z, in eq. (|16|) vanishes and i^^' 
becomes equal to -F,r- Thus in eq. H33() we have replaced F^/ with i<V and at the same time Q 
with Q, although in this limit the integral of the divergence of the singlet axial current does 
not vanish. Recalling the expression (|12j) and the index theorem (jl3|) . one can equivalently 
write for the (quenched) rj' mass the suggestive formula 



n 



2Nf 

where V is the space-time volume of the lattice 



= lim ^^"-^ - "-^^^ . (34) 
Nf/No=0 V^oo V 



^For a discussion of several subtleties on the derivation of this formula see ref. 
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3. Paying due care to flavour matrix normalization, one can combine the non-singlet WTI (writ- 
ten for a given fixed flavour, h) 



= 2m / fx {P''{x)P^{0)) - -^(5°(0)) (35) 

Nf 



with eq. obtaining 



/f zv 
d^x (Q(x)Q(O)) = (2m)2 J d^x (pO(x)pO(O)) | , (36) 

where the superscript ZV means that only Zweig- Violating (hairpin) diagrams should be 
included in the r.h.s. of the equation in carrying out the fermion functional integral. For GW 
fermions the relation H36|) is an algebraic identity that can be directly proved, using eq. @, 
the definition (|12|) and the explicit expressions of P° and S° (eqs. ^ and A few 

observations are in order here. 

• An equation like ()36() holds also if the fermion determinant is neglected (quenching) . This 
follows immediately from to the GW relation, after the fermion integration is performed, 
as shown in the Appendix. 

• In the quenched limit the r.h.s. of eq. (|36() possesses a double pole at vanishing quark 
(pion) mass with a residue related to m^, JSIEI- To be precise its residue is in our 
normalization m^;m^ F^/2Nf This observat ion was the basis of the many quenched 
simulations carried out in the years in lattice QCD aimed at extracting the mass of the 
pseudo-scalar flavour singlet, starting from the seminal work of Hamber, Marinari, Parisi 
and Rebbi (see refs. [13 ED 1221) • 

• With the idea of trying to set up a formula for the r]' mass which would not depend on 
the details of the lattice definition of the topological charge density operator, in ref. 
for Wilson and few years later for staggered fermions [21], the equation 



2Nf 



ZV 

(37) 

quenched 



was argued to hold in the limit Nj/Nc — > ^. Eq. (|37() can be regarded as an expression 
of the residue of the double pole present in the / d'^x {P^{x)P^{0)) correlator at 

Nf/Nc = (i.e. in the absence of the fermion determinant). In this way no use of the 
identity (|36|) is actually made, though eq. (|37|) is obviously consistent with (|36]) for GW 
fermions. 

4. As in the formal continuum theory ^S], the chiral condensate can be extracted from the small 
m expansion of the (lattice) full topological susceptibility, as defined by eqs. (|12|) and (|22]) 



4iV/xIf = -4m(50(0)) +0(^2), (38) 

•' m=U 



*We are using the definition {0|9fi.451 1?;') = y/2NfF^'m'p with the identification of F^/ with in the quenched 
hmit. We are also assuming that the pion bound state exists even in the quenched theory. 

^ActuaUy the way the r;' mass formula was written in ref. |2.-i) is wrong by the finite normalization factor 1 — 
dM°{M)/dM\M^M,r, see ref. ^ for notations. 
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with no need of power divergent subtractions. 

5 Conclusions 

In this paper we have shown that there are no m? /a^ power divergences in Xjl') defined in 
eq. (P^. Thus, if GW fermions are used, the topological susceptibility of the full theory 



can be employed to extract the physical value of the chiral condensate using the eq. , without 
the need of performing any power subtraction. The reason is that the power divergence in {S^), 
outside the chiral limit, is exactly compensated by a similar divergence in the (P^P^) correlator in 
the r.h.s. of eq. (|22|) . 

Actually we have proved more than that. We have proved that all power divergences in the 
(P^P^) correlator arise from the ZC contributions. This conclusion follows from the non-singlet 
WTI (|35|) . where we see that only ZC diagrams contribute. 

This observation may be of some interest in view of the formula (jSZI), where the (quenched) 
r]' mass is expressed in terms of the ZV quenched correlator of two singlet pseudo-scalar quark 
densities. In this case, in order to conclude that power divergences are also absent from the r.h.s. 
of eq. (|S7|). one must assume that taking the limit Nf/Nc ^ in eq. is equivalent to dropping 
the fermion determinant and does not introduce unexpected 1/a^-power divergences. If this is the 
case, one can imagine to check the validity of the formula ()33() and the assumptions underlying it by 
comparing the value of the YM topological susceptibility (r.h.s. of eq. (jSU) with the residue of the 
double pole that arises at zero quark mass when the fermion determinant is dropped from eq. ()36() 
(thus ending up with eq. (|37j) ). In this context it is interesting to mention that the YM topological 
susceptibility has been recently computed at several values of the lattice spacing by counting the 
number of zero modes of the Neuberger-Dirac operator [27\ I28j . Data are compatible with the 
scaling behaviour expected for a quantity of dimension d = A and no sign of power divergences 
(within errors). 

Appendix 

In this appendix we want to prove that the formula 



where ((. . .)) means fermion field contraction only, holds configuration by configuration in the 
massive theory. To this end we recall that the expression of the massive GW operator in eq. © is 




(39) 




(40) 



Dra = {l-^)D + m. 



(41) 
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In terms of Dm the GW relation @ takes the form 

-I^Dm + Dmlh - 2m75 = _"amx ij^mlbDm " rn{'-I^Dra + Dmlb) + m^lb) ■ (42) 

\^ 2 I 

Putting similar terms together and multiplying by the inverse of the massive GW operatore, D^, 
we get from ()42() . after taking the trace and integrating over space-time 

-mj d^xTr b^D-^] = ^J d^xTv [-/^Dm] • (43) 

At this point it is enough to observe that 

1) contracting the fermion fields in the operator appearing in eq. (|4()j) gives 



2) tracing Dm with 75 yields 



(44) 



nrr) 

TT[j,Dm] = {l-—)Tt[^,D]. (45) 
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